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Quantum criticality in the 
two-channel pseudogap Anderson 
model: A test of the non-crossing 
approximation 
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We investigate the dynamical properties of the two- 
channel Anderson model using the noncrossing 
approximation (NCA) supplemented by numerical 
renormalization-group calculations. We provide ev- 
idence supporting the conventional wisdom that the 
NCA gives reliable results for the standard two-channel 
Anderson model of a magnetic impurity in a metal. We 
extend the analysis to the pseudogap two-channel model 
describing a semi-metallic host with a density of states 
that vanishes in power-law fashion at the Fermi energy. 



1 Introduction Quantum criticality is currently pur- 
sued across many areas of correlated matter, from insu- 
lating or weak magnets to metal-insulator systems to un- 
conventional superconductors. Interest in continuous phase 
transitions at temperature T = is motivated both by the 
appearance of novel phases near such transitions and by 
a richness of quantum critical states that extends beyond 
the traditional description in terms of order-parameter fluc- 
tuations. A fertile area for experimental study has been 
the border of magnetism in rare-earth and actinide inter- 
metallics, where there is mounting evidence that a cru- 
cial issue is the fate of the Kondo effect on approach to 
the quantum phase transition HE): if Kondo screening be- 
comes critical concomitantly with the bulk magnetization, 
then dynamical scaling ensues, signaling the absence of ad- 
ditional energy scales at criticality J3]|4]. 



This model exhibits continuous quantum phase transi- 
tions between weak- and strong-coupling phases. The 
NCA is shown to reproduce the correct qualitative fea- 
tures of the pseudogap model, including the phase dia- 
gram, and to yield critical exponents in excellent agree- 
ment with the NRG and exact results. The forms of the 
dynamical magnetic susceptibility and impurity Green's 
function at the fixed points are suggestive of frequency- 
over-temperature scaling. 
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This paper investigates critical Kondo destruction and 
dynamical scaling in the pseudogap two-channel Anderson 
model. We focus on the treatment of the problem using 
the noncrossing approximation (NCA), which involves a 
self-consistent evaluation of all irreducible self-energy di- 
agrams without vertex corrections. Within a pseudoparticle 
representation of the impurity spin [5 1, the NCA threshold 
exponents of the T = pseudoparticle propagators for the 
one-channel Anderson model J6) are known to differ from 
the correct exponents inferred from Friedel's sum rule [f7]. 
In contrast, the NCA threshold exponents for multichannel 
Anderson models agree |8] with those predicted by bound- 
ary conformal field theory [9.jl0|. The NCA is therefore 
believed to work well for this model, which has been stud- 
ied extensively as a relatively simple route to non-Fermi 
liquid behavior. The pseudogap version of the two-channel 
Kondo model, in which the host density of states vanishes 
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in power-law fashion at the Fermi energy, was first stud- 
ied in Ref. ifTTIl . The model and the counterpart Anderson 
model have attracted renewed interest [:l2| due to propos- 
als Ml 311 141 of their realization in the context of adatoms 
on graphene, where the band structure gives rise to two 
symmetry-inequivalent Dirac points. Recent scanning tun- 
neling spectroscopy data are in line with expectations for 
pseudogap multichannel Kondo physics fl5l . 

Here we briefly revisit the reliability of the NCA for 
the metallic two-channel Anderson model before turning 
to the pseudogap case, where we apply a scaling ansatz to 
obtain T = threshold exponents, and hence extract criti- 
cal exponents describing physical properties. We find good 
agreement between these exponents and ones obtained us- 
ing full numerical solutions of the NCA equations at T > 
and using the numerical renormalization group. 

2 The pseudogap M-channel Anderson model 

The SU(iV) xSU(M) Anderson model can be written 

N M 

hmca = X X e ^ c L„ c ka„ + e f X fi>f° 

k <j— 1 /i= 1 cr 

+ V J2 (i) 

where creates a conduction electron of wave vec- 
tor k in channel /i and spin projection a. The repre- 
sentation dj. = of the impurity electron creation 
operator in terms of pseudofermion creation combined 
with slave-boson annihilation faithfully reproduces the 
SU(7V)xSU(Af) Anderson model provided that the con- 
straint Q = J2a ftfa + E M b t b t, = 1 is enforced exactly. 
We assume a conduction-electron density of states 

p( U ) = -±TmG c (u) = £^\w\ r e(D-\w\) (2) 

that vanishes at the Fermi energy u> = in a manner gov- 
erned by exponent r, taken to satisfy < r < 1; the spe- 
cial case r = describes a flat (metallic) band. The half- 
bandwidth D acts as the basic energy scale in the problem. 

N = 2 and M = 1, Eq. (fl3 reduces to the standard 
Anderson impurity model (with infinite on-site Coulomb 
repulsion) which in turn contains the spin-isotropic Kondo 
model as an effective low-energy model, while for N — > oo 
various saddle-point approximations can be constructed 
|[T6l[8lfT7ll . The dynamical large-iV limit of Ref. E) for the 
spin-isotropic Kondo model uses a generating functional 
equivalent to that for the NCA, making it seem natural to 
expect similar results from the two approaches. However, 
it is important to note that within NCA, the slave-boson 
propagator Gb(w) is not a Hubbard-Stratonovich decou- 
pling field and is therefore dynamic even at the bare level, 
and that Gb(uj) couples to the local constraint, which is 
enforced exactly. As a consequence, Gbioj) and the pseud- 
ofermion propagator Gf(uj) in the NCA develop thresh- 
old behavior reminiscent of the core hole propagator in the 



x-ray edge problem, resulting in maximally particle-hole 
asymmetric spectral functions: ImGf,(w — > 0,T = 0) ~ 
0(uj)\uj\- ab and Im G/(u; -> 0,T = 0) ~ G(ui)\u)\- a ? . 
Therefore, G b (r, T) and G/(r, T) should be rather differ- 
ent from their dynamical large- N counterparts, raising the 
question of whether dynamical or ui/T scaling — a prop- 
erty that arises naturally within the dynamical large-iV 
approach — can carry over to the NCA. 

3 Asymptotically exact zero-temperature so- 
lution For the case r = of a constant density of 
conduction-electron states, an exact finite-temperature 
solution can be obtained by transforming the NCA's in- 
tegral equations into a set of coupled differential equations 
J6j. As this procedure relies on specific properties of the 
NCA solution for r = 0, its extension to the pseudogap 
case is unclear. We can address the zero-temperature form 
of the NCA solution by imposing a scaling ansatz for 
the pseudoparticle spectral functions. The scaling ansatz 
has been successfully used to extract critical properties 
of impurity models treated within the dynamical large-iV 
approximation, including the Kondo model both with a 
metallic (r = 0) density of states IfTTl and with a pseu- 
dogap flin . the Bose-Fermi Kondo model fl9l , and the 
pseudogap Bose-Fermi Kondo model l20l . Here we show 
that a generalization of this ansatz to the case of extreme 
particle-hole asymmetry (generated by the exact enforce- 
ment of the constraint) can be used to extract the critical 
properties of the fixed points within the NCA. 

The NCA self-energies for the pseudoparticle propaga- 
tors at real frequencies are H2 1 11221 

2%(u>) = V 2 Y, I ^W(-e) G?*(e+w), (3) 

ATM = y2 E / ^W £ ) G / ot ( £ + w )> ( 4 ) 

where the superscript "ret" specifies a retarded function, 
A cail (u) = -7r- 1 ImG^(o;) = p(ui), and /(e) is the 
Fermi-Dirac distribution function. Working at temperature 
T = 0, we make the ansatz 11231171 

A p (oj) = a p <9(cj)w~ Qp , p = f,b, (5) 

for the pseudoparticle spectral functions A p (oj) = —tt^ 1 
lmG p ct (u)), with a p being a threshold exponent and a p a 
constant. Substituting these scaling forms into Eqs. (f3]i and 
(|4), and applying Dyson's equations, one obtains the self- 
consistency conditions 

aj + a b = 1 + r, (6) 
sin (jraf) sin [tt (r — a/)] x 
f M sm[n (r - a f )} | sin( 7 ra / ) \ Q (7) 
\ N ctf 1 + r — ctf ) 

Equation (O is obtained by matching exponents and Eq. 
(UJ by matching amplitudes. Hereafter, we focus on the 
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Figure 1 Left: Variation with pseudogap exponent r of 
possible solutions a/ of Eq. (0 for the two-channel An- 
derson model. Dashed lines represent unstable fixed points 
and dotted lines show values not observed in numerical 
solutions of the finite-temperature NCA equations. Right: 
Renormalization-group (RG) flows and fixed points as 
functions of hybridization V for the ranges < r < ro 
and r > ro. Labels (I)-(V) connect RG fixed points with 
NCA solutions in the left panel. 
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Figure 2 Pseudoparticle spectral functions Af(u>) and 
At (uj) for a metallic host (r = 0), calculated at tempera- 
ture T = 5xlO- 8 Dfore f /D = -0.6 and (V/D) 2 = 0.5. 
The NCA pseudoparticle threshold exponents otf = ctb = 
0.5 agree with the exactly known results. 



two-channel case N = M = 2, for which the possible so- 
lutions of Eqs. © and © are plotted schematically on the 
r-ctf plane in Fig. [I] For r = 0, the solutions are aj = 0, 
otf — 1 (local moment), and Q-f = h (intermediate cou- 
pling), in agreement with Ref. [7|. For < r < ro, where 
the condition (ro + l)7r/2 = cot(ro7r/2) yields anumerical 
value ro ~ 0.292, there are five solutions, of which three 
correspond to stable renormalization-group fixed points: 
local moment (I), two-channel Kondo (III), and infinite-c7 
resonant level (V). For r > ro, there are just three solu- 
tions, of which only (I) and (V) are stable. 

The (physical) impurity Green's function is obtained as 

G£» = J dee-^[G%(e+u;)A bfl (e) 

-Gff(e- U )A f r(ej\, (8) 

where f3 = 1 /T and "adv" means advanced. The T = 
scaling form of this Green's function and of the local spin 
susceptibility \ can be deduced from the pseudoparticle 
propagators by simple arguments, yielding ImGj ct (u) oc 
\uj\ 1 ' a f~ a " ee \uj\~ r andlmx(w) oc sgn(w) l^l 1 " 2 ^. 

The NRG has been applied previously to the pseudogap 
two-channel Kondo model 11 10121 . For r < r max ~ 0.23 
(somewhat smaller than the NCA value ro ~ 0.292), the 
method finds three stable fixed points separated by two un- 
stable critical points, precisely analogous to the situation 
shown in Fig. Q] For r max < r < 1, the NRG yields one 
critical point between two stable fixed points, again in one- 
to-one agreement with the NCA scaling ansatz. Of these 
fixed points, only those corresponding to (I) — (III) can be 
accessed for level energies ef < 0, the regime considered 
in the present paper. (The NCA treatment of cases with 
ef > will be reported elsewhere 11201 .) 

We have confirmed that the pseudogap two-channel 
Anderson model shares the same fixed points and criti- 
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Figure 3 Static susceptibility \ vs temperature T for a 
metallic host (r = 0) and £//£>= -0.6, (V/D) 2 = 0.5. 
The NCA captures the logarithmic divergence of x(T)- 



cal properties as its Kondo counterpart, in agreement with 
the conclusions of Ref. [12|. This justifies our compar- 
isons below between NCA results for the Anderson model 
and NRG results for the Kondo model (the smaller Hilbert 
space of which allows greater numerical efficiency). 

4 Finite-temperature NCA solution At tempera- 
tures T > 0, the NCA equations are amenable to a numer- 
ical solution on the real frequency axis over a wide param- 
eter range. Details of the numerical evaluation scheme can 
be found, e.g., in 1241 . 

Results for r = 0: As mentioned above, the NCA pre- 
dicts the correct threshold exponents for multichannel An- 
derson models with a metallic density of states. This is il- 
lustrated in Fig. |2] where the fitted pseudoparticle expo- 
nents agree very well with the NCA scaling ansatz and with 
the boundary conformal field theory for the two-channel 
Kondo model (25! . Figure [3] demonstrates that the NCA 
correctly captures the logarithmic divergence in temper- 
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(I) 


(II) 


(III) 


(IV) 


(V) 


Of 


1.0 


0.967 


0.575 


0.183 


0.15 


a b 


0.15 


0.183 


0.575 


0.967 


1.0 



Table 1 Pseudoparticle threshold exponents at T = for 
the two-channel Anderson model with r = 0.15. 




Figure 4 Pseudoparticle spectral functions Af{uS) and 
A b (u) for r = 0.15, T = 5 x 10~ 8 D, e f /D = -0.6, 
and (V/D) 2 = 0.27. The power-law behaviors correspond 
to solution (II) in Table Q] and Fig.Q] 




Figure 6 Impurity spectral function A d (u)) for r = 0.15 
very close to the critical fixed point (II) in Fig. Q] calcu- 
lated using the NCA at temperature T = 5 x 10~ 8 .D for 
e//D = -0.6 and (V/D) 2 = 0.27. Also shown is the 
spectral function calculated within the NRG for T = 0, 
e f /D = -0.1 and (V/D) 2 ~ 0.0486. The NCA and NRG 
spectral functions are described by very similar exponents 
(0.166 and 0.150, respectively), and both show that at the 
critical point, A d (uj) is particle-hole symmetric at low en- 
ergies. 




co/T 



Figure 5 Dynamical scaling of x(w, T) at the critical fixed 
point (II) in Fig. OQ for r = 0.15, e f /D = -0.6, and 
(V/D) 2 = 0.27. 

ature T of the local spin susceptibility \- The sublead- 
ing behavior of the impurity spectral function, A d (u) = 
-ir^ImGf, is also reproduced 0: A d (ui) - A d (0) ~ 
\J\u>\. Taken together, these pieces of evidence indicate 
that the NCA gives qualitatively correct results for the 
r = two-channel Anderson problem. 

Results for r = 0.15: We now turn to numerical results 
for the two-channel pseudogap Anderson model, focusing 
first on the case r = 0.15 with e/ = — Q.6D as a represen- 
tative example of the behavior in the range of pseudogap 
exponents < r < Tq. For this case, Eqs. © and (0 
predict five solutions, listed in Table Q] Solution (I) cor- 
responds to the local-moment fixed point where V effec- 



tively vanishes. Obtaining converged numerical solutions 
near this fixed point is very difficult as the resulting sharp 
features cannot be resolved on discrete frequency grids. In- 
creasing the hybridization V until a solution can be stabi- 
lized at the lowest accessible temperatures yields the crit- 
ical solution (II) at V = V c where (V c /D) 2 ~ 0.27. Fig. 
|4] shows that the pseudoparticle exponents at this critical 
point agree well with the scaling ansatz results in Table Q] 
The temperature dependence of the static susceptibility re- 
flects the frequency behavior. As a result, the dynamical lo- 
cal spin susceptibility T) at the critical point exhibits 
the dynamical scaling form (see Fig. |5j 

X (w,T) = T-*$(w/T), (9) 

with x = 0.934. It is instructive to compare this result 
against other methods. The NRG is unable to reliably ac- 
cess the regime < M/7 1 <C 1, so it cannot fully test for 
dynamical scaling. However, we find for the two-channel 
Kondo model with r = 0.15 that at the critical point, 
the NRG gives x(w = 0,T) cx T~ x and Imx(w,T = 
0) cx \u\~ y with x = y = 0.930 ± 0.001. These proper- 
ties are entirely consistent with Eq. (0, and show that the 
NCA does an excellent job of calculating the exponent x. 
We note that this exponent deviates significantly from the 
leading-order value x = 1 — 2r 2 = 0.955 coming from an 
expansion about the local-moment fixed point Ifl2l . 

Figure [6] shows the impurity spectral function A d (uj) 
very close to the critical fixed point (II) in Table Q] The 
power-law variation is compatible with the |w|~ r expected 
from the scaling ansatz, which coincides with the exact 
behavior known to hold at all intermediate-coupling fixed 
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Figure 7 Pseudoparticle spectral functions Af(oj, T) and A b (uj, T) at the fixed point (III) in Fig.Q]for r = 0.15, e//D = 
—0.6, and (V/D) 2 = 2.0. (a) Frequency variation at T = 5 x 10 -8 . (b), (c) Dynamical scaling at different temperatures. 
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Figure 8 Dynamical scaling of the impurity spectral func- 
tion Ad(io, T) at the stable fixed point (III) in Fig. Q]for 
r = 0.15, tf/D = -0.6, and (V/D) 2 = 2.0. 




Figure 9 Static susceptibility \ vs temperature T at the 
stable fixed point (III) in Fig. [TJ for r = 0.15, ef/D = 
-0.6, and (V/D) 2 = 2.0. 



points ifPH and with the NRG results also plotted in Fig. 
Both the NCA and the NRG show the critical spectral 
function to be particle-hole symmetric at small \uj\. 

The stable fixed point (III) in Table [TJ can be accessed 
by increasing the hybridization beyond the critical value 
V c . Figure Qpresents the pseudoparticle spectral functions 
at various temperatures for (V/D) 2 = 2.0. The threshold 



exponents a/ = 0.58 and ctb = 0.575 extracted from Fig. 
ITJa) are in line with the prediction af = ctb = (1 + r)/2 of 
the scaling ansatz. Figures |7Jb) and [3c) demonstrate that 
the temperature dependences of Af(oj,T) and Ab(u>,T) 
are compatible with the frequency behavior such that dy- 
namical, or w/T-scaling ensues. This carries over to the 
impurity spectral function, which, as seen in Fig. [8] is com- 
patible with the scaling 



GT(u,T)=T- r V(uj/T) 1 



(10) 



This scaling is consistent with the exact result lfl2l men- 
tioned above, i.e., Ad(oj, T — 0) oc M~ r . 

Figure [9] plots the static local susceptibility vs temper- 
ature at the stable fixed point (III). A rather narrow win- 
dow of asymptotic temperature dependence — presumably 
a consequence of a strong subleading contribution to 
l/x(w = 0,T) — gives an exponent x = 0.195 in rea- 
sonable agreement with the scaling ansatz prediction 
x = 2a f - 1 = 0.15. 

Results for r = 0.4: We end this section by consider- 
ing a representative case in the range rrj < r < 1. For 
r = 0.4, the scaling ansatz predicts an asymmetric crit- 
ical point described by ay = at, = (1 + r)/2 = 0.7, 
the solution corresponding to (III) in Fig. Q] Figures [TOl 
and [9] show results in the vicinity of this critical point, ob- 
tained for e f = -0.55D and (V/D) 2 = 5.8. The pseu- 
doparticle exponents a/ = 0.7, a& = 0.68 extracted from 
Fig. [Tol a) closely follow the scaling ansatz, and the dy- 
namical scaling of Ad(uJ,T) with an exponent 0.37 ~ r 
[Fig. [T0l b)l agrees with the exact result JTZj. Furthermore, 
the exponent y = 0.44 of x(w,T = 0) oc \lo\ y fitted 
from the rather narrow frequency window of power-law 
behavior in Fig. [TOl b) is in reasonable agreement with the 
value y ~ x ~ 2af — 1 = 0.4 predicted by the scaling 
ansatz under the assumption of dynamical scaling and with 
x = y = 0.381 ± 0.001 given by the NRG. The impurity 
spectral function, shown in Fig. [TOl c). also appears to be 
consistent with dynamical scaling. 

5 Conclusion We have investigated the reliability 
of the non-crossing approximation (NCA) for the two- 
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Figure 10 Results for r = 0.4 in the vicinity of the asymmetric critical point (III) in Fig.Q] calculated for ej — — 0.55D 
and {V/D) 2 = 5.8. (a) Pseudoparticle spectral functions Af(u,T) and A b {u,T) at T = 10~ 8 D. (b) Dynamic spin 
susceptibility xi^i T) at T = 10~ 8 D. (c) Dynamic scaling of the local density of states Ad(w, T). 



channel pseudogap Anderson model. This was accom- 
plished by comparing finite-temperature, finite-frequency 
solutions of the NCA equations with asymptotically ex- 
act zero-temperature NCA solutions, with numerical 
renormalization-group calculations, and with exact results 
where available. In contrast to the well-known shortcom- 
ings of the NCA for the single-channel Anderson model 
with a constant density of states at the Fermi energy, the 
NCA captures surprisingly well the asymptotic low-energy 
properties of the two-channel model, both for metallic and 
semi-metallic (pseudogapped) hosts. In cases of a pseudo- 
gap, the results that we have presented for the magnetic 
susceptibility and the impurity spectral function are sug- 
gestive of frequency-over-temperature scaling in the dy- 
namical properties. More complete testing for dynamical 
scaling is planned as future work. Finally, we note that the 
validation of the NCA treatment of this problem at equi- 
librium opens the way for its extension to nonequilibrium 
steady-state conditions l20l . 
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